Detecting non-Abelian geometric phase in circuit QED 
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We propose a scheme for detecting noncommutative feature of the non-Abelian geometric phase in circuit 
QED, which involves three transition qubits capacitively coupled to an one-dimensional transmission line res- 
onator. By controlling the external magnetic flux of the transmon qubits, we can obtain an effective tripod 
interaction of our circuit QED setup. The noncommutative feature of the non-Abelian geometric phase is man- 
ifested that for an initial state undergo two specific loops in different order will result in different final states. 
Our numerical calculations show that this difference can be unambiguously detected in the proposed system. 
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Phase factor has played a profound role in quantum physics. 
Apart from the familiar dynamical phase, geometric phase 
(GP), discovered in 1984 by Berry JT|, has deep physical 
meanings. Berry pointed out that GP occurs when a sys- 
tem is subjected to a cyclic adiabatic evolution, which results 
from the geometrical properties of the parameter space of the 
Hamiltonian. Especially, Wilczek and Zee found non-Abelian 
gauge phase results from Berry's formula in 1984 |0. GP de- 
pends only on the solid angle enclosed by the parameter path, 
and thus is robust against local noises Jirjgt]. Therefore, it 
has been proposed to implement fault-tolerant quantum logi- 
cal gates for universal quantum computation Ill0l4l5ll . 

Non-Abelian GP differs from Abelian geometric phase ac- 
cording to the commutation features of their gauge potential, 
i.e., an initial state undergoes two specific cyclical evolution 
in different order will result in different final states. Up to 
now, AbeljanGP has been experimentally detected in various 
systems Jlrjl- 181. while non-Abelian GP has not been verified 
yet. Usually, dark states are used in detecting GP so that dy- 
namical phase will not appear. Conventionally, it is the tripod 
Hamiltonian that has been proposed to detect the non-Abelian 
GP M19H21I1 . However, it is difficult to find tripod configura- 
tion atomic energy levels, which impedes experimental detec- 
tion of the non-Abelian GP. Recently, it is also proposed that 
non-Abelian GP can be detected by two laser beams interact- 
ing with a three-level A atom in cold atomic system B221I23I1 . 
Since two of the eigenstates in this scheme are only near de- 
generate, dynamical phases will also occur during the pro- 
cess, and thus one needs additional effort to conceal it 12211 . 
Meanwhile, the non-Abelian GP is also proposed to be de- 
tected in a new designed multi-level superconducting circuit 
02411 . However, multi-level scenario of this superconducting 
nanocircuit is very sensitive to its background charge noise. 
Therefore, to certify the fundamental non-Abelian nature of 
the non-Abelian GP, it is of great importance to find an ex- 



perimentally accessible system that can host the exotic non- 
Abelian structure. 

Superconducting system is regarded as one of the most 
promising candidates for physical implementation of qubits 
which can support scalable quantum information processing 
II25[ 12611 . Furthermore, by placing su perc onducting qubits 
in a cavity, i.e, circuit QED setup J27, 28], the system will 
have several practical advantages including strong coupling 
strength, immunity to noises, and suppression of spontaneous 
emission. Here, we propose to detect the noncommutative fea- 
ture of non-Abelian GP with effective tripod Hamiltonian in 
circuit QED. The setup we consider consists three transmon 
qubits that are capacitively coupled to an one-dimensional 
(ID) high-Q transmission line resonator (cavity), which has 
recently been realized experimentally n29ll . With proper cho- 
sen parameters, such setup can be effectively described by the 
tripod Hamiltonian 13011 . and thus can be used to detect the 
noncommutative feature of non-Abelian GP. Furthermore, the 
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transmon qubit possesses remarkable superiority [3JJ, e.g., it 
achieves exponential insensitivity to charge noise without in- 
creasing the sensitivity to either flux or critical-current noise. 
Note that when adding a shunt capacitor to a flux qubit will 
also lead to low-decoherence qubit l32ll . 

The considered circuit QED architecture is shown in Fig. [T] 
with three identical transmon qubits that are capacitively cou- 
pled to the cavity. The transmon qubit has effective Josephson 
energy Ej = Ej tmax \cos(ir$/$ )\ with Ej tmax , $ and $ 
being the Josephson energy of the Josephson junctions, the ex- 
ternal magnetic flux and the flux quantum, respectively. This 
type of qubit has good coherence performance. The charg- 
ing energy of the transmon is much small compared with the 
Josephson energy {E c <^E.j). With Ej/E c = 50 (E c = 0.3 
GHz, Ej = 15 GHz), the energy difference of the two low- 
est levels (defined as first excited state |e) and ground state 
\g)) is approximately ^J&EjEc, and the relaxation time for 
|e) is on the order of 0.06 s Bill . For an ID cavity with length 
L = A = 1 cm, we can get rms voltage V® ms — y hui/lc of an 
antinode between two superconducting lines, where I and c are 
the inductance and capacitance per unit length, respectively. 



2 





Z = A = l cm 




\ 













1/ 
















FIG. 1: (Color online) Schematic of the circuit QED architecture. 
The ID superconducting line resonator consists of a full-wave sec- 
tion ( L = A = 1cm ) of superconducting coplanar waveguide. 
Qubits placed between superconducting line resonators consist of 
two small Josephson junctions, which permit tuning of the effective 
Josephson energy by an external flux $ in a 1 /im loop. There is 
approximate 10 /im gap between superconducting line resonators. 
The input and output signals can be coupled to the resonator via the 
capacitive gaps in the center line. 



|0) and |1) denote the cavity mode has and 1 microwave 
photon. Then, in this subspace, the interaction Hamiltonian 
([TJ can be written as 
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(3) 



where Aj = e< — oj is the detuning of ith qubit of the cav- 
ity. Here, we choose to drive the system by means of a 
small time-dependent quantities and separate them from the 
time-independent system denoting with superscript (0). As 
gi and Ai are both related to <£>i, the time-dependent driven 
can be added by choosing the magnetic flux as $j(t) = 
< f>- ' 1 + <5$i(t), and the corresponding Hamiltonian can be 
written as H su b — H^> + 8H(t). Assuming that the time- 
dependent fluxes 6$>i(t) oscillate with the frequencies oji/2'K, 
the corresponding quantities are written as 



d®i(t) = Fi(t) cos[ojit + cpi 



(4a) 



As a result, qubits are coupled to the superconducting line by 
means of the voltage V = V® ms (a + a'). Remarkably, for 
coplanar waveguide cavity, cavity quality factor Q ~ 10 6 has 
already been demonstrated ll33ll . which means that the inter- 
nal losses can be very low. With three qubits fabricated at the 
antinodes of the cavity voltage, the strength of the coupling 
to the resonator is maximized for all three transmom qubits. 
Then, the system can be described by the Tavis-Cummings 
Hamiltonian 
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where uj = \j\fTc is the resonator frequency, a (aJ) 
is annihilation (creation) operator of the ID cavity mode, 
Si = y/8E£Ej/h is the energy splitting of the ith qubit, 
a 1 are Pauli operator for the ith transmon, and = 
2eC n g V° ms (E l J /8E i c ) 1 / i /V2Ci is the strength of coupling 
between the ith qubit and the superconducting line. For 
C g /C s = 0.1 with C E = C g + 2Cj, the coupling strength 
controlled by external magnetic flux <£>i will be on the order of 
100 MHz M28I1 . Describe e% and gi of Hamiltonian in Eq. (UJ 
in external magnetic flux they read 



M*i) = \/8J^.mJ«»( 
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Explicitly, the coupling strength gi and the qubit frequency e, 
are endowed with a relation gi($>i) oc y/ei($i). 

Restrict the system in the only one-excitation subspace of 
the dynamics states {\lggg), \0egg), \0geg), \0gge)}, where 



6Ai(t) = Li(t) cos[uJit + tpi], 



Sgi(t) = Ti(t) cos[cjjt + <pi 



(4b) 



(4c) 



where amplitudes Li (t) and Tj (t) are determined by means of 
externally modulated flux amplitudes Fi (t) based on equation 
55 and d2bl 

The eigenvalues of the main Hamiltonian in the one 



excitation subspace are {0, A]/ , A 



(0) A (0) 



A 3 0) }. 



In this eigen- 



basis, choosing 8gi(t) and 5Ai(t) to oscillate with frequency 
Ui = A\° , the effective Hamiltonian in rotating frame with 



rotating wave approximation reads [30] 
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where effective Rabi frequencies are 



(6) 



with time-independent parameter r\i 

ff f)/(2Af). 

To detect the noncommutative feature of the non-Abelian 
GP, we now parameterize Rabi frequencies in Eq. (|6]l to form 
two specific evolution loops C\ and C2 with U\ and U2 being 
their respective evolution operators. The non-Abelian nature 
of GP is verified by the fact that for an initial state undergoing 
the two specific loops in different order will result in different 
final states. This noncommutative feature of the gauge struc- 
ture leads to the non-Abelian characteristic of the non-Abelian 
GP. For convenience, the initial phase of Rabi frequencies are 
chosen as (fii — if2 = 0, ^3 = £, respectively. Therefore, 
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only the amplitudes of Rabi frequencies vary with time, which 
are determined by the amplitudes -Lj(t) of the external mag- 
netic flux. Modulate Li (t) appropriately so that the two loops 
C\ and C 2 are obtained as 



d : n x = n /(t), 
c 2 : ni - n f(t), 



n' 2 = an f(t), 



n' 3 = floe 
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where f(t) = cos[7rf/(2r)] for an interval of t E [— r, r] and 
£ = irb/r. Two variables a and /3 make a distinction between 
the loops C\ and C 2 with /? being a time delay factor. At prac- 
tical parametrization, Q 2 = f2 / 2 (i) in the loop C\. To form 
fl' 2 in loop C 2 , we introduce another magnetic flux which will 
produce Vl± = (a — 1)Hq f 2 (t), which is turned on only when 
forming the loop C 2 . 

Then, Rabi frequencies f2j in C\ can be rewritten as 



Sli = il sin cos ^e^ 1 . 
fi 2 = O sin sin 0e^ 2 , 
Q 3 = flcosOe^ 3 



(8) 
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Then, two dark eigenstates of H in Eq. (0) are 
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where {|1), |2), |3)} denote {\0egg) : \0geg), \0gge)}, respec- 
tively. We then can get the gauge potential A based on 
.1,,,, I), + I), as if 
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Therefore, A\ is 
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and its corresponding time evolution operator is 



E/i=Pexp(- / Ai^dx^), 
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FIG. 2: (Color online) (a) Schematic of the population difference Pd 
varying with parameters a and ft, (b) The populations P (red), P' 
(black) and Pd (blue) as functions of the parameter (3 (0 — ^ 1.0) for 
the parameter a = 6. 



where V denotes the path-order operator. In order to unam- 
biguously detect non-Abelian geometric phase, we confine pa- 
rameters (8, <fi, £) vary from (0, 0, — tt) to (0, 0, it) with time 
— r — > t. Similar to loop C\, we can get the corresponding 
evolution operator U 2 based on Eq. dTbl in the loop C 2 . 

To detect the non-Abelian nature, we first prepare the ini- 
tial state as \ip)i = \D 2 )i = |1), and then let it undergo two 
closed paths C\ and C 2 in different orders, i.e., U = U 2 Ui 
(first Ci then C 2 ) and U' = U X U 2 (first C 2 then d). In or- 
der to implement the evolution U(U'), let XJ\ (U 2 ) in effect 
during time — r — > r, while f/ 2 (L^) during time r — > 3r. 
The final states will be \ip) f = ?7|1) = -J7 2 i|2) + £/ 22 |l), 
|^ = U'\l) = -Z7 21 |2) + ?7 22 |1), respectively. Therefore, 
the population difference Pd of the two different final states in 

ID is 



P d = P' -P = \U' 



\u 2 



(14) 



Whenever the Pd 7^ is detected, the noncommutative feature 
of the non-Abelian GP is verified. The population difference 
Pd is numerically calculated with variables a and /3, as shown 
in Fig. |2](a), which obviously indicates that Pd 7^ 0. Fig. |2](b) 
is a specific plot of the population difference Pd with j3 as the 
only variable while a = 6, which shows maximum Pd ~ 0.6 
when /3 = 0.9. 
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Detecting the population difference of state |1) = \0egg) 
means that we just need to observe population difference on 
the excited state \e) of qubit 1, which can be realized by 
quantum non-demolition (QND) measurement. This can be 
achieved by tuning the qubit dispersively coupled to the cav- 
ity with a large detuning A, and then measuring Hamiltonian 
will be Hm = h(u> + x <J z)a !< a with \ = <7 2 / A. We can then 
get a different frequency shift ±x of the cavity mode with 
the qubit state on \g) and |e), respectively. With the coupling 
strength g\/2-K = 100 MHz and the detuning A = 5g, we can 
get the frequency shift as \/^ n = 20 MHz, which is readily 
resolvable experimentally with high fidelity [|28ll . 

In summary, we have proposed an experimentally feasi- 
ble scheme to detect the noncommutative feature of the non- 
Abelian GP with effective tripod Hamiltonian in circuit QED. 
The non-Abelian nature of GP is verified by the fact that for an 
initial state undergoes the two specific loops in different order 
will lead to different final states. This differences is detected 
through observing the population difference of state \0egg), 
which is achieved QND measurement in circuit QED. 
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